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We discuss different approaches for studying the influence of disorder in the 
three-dimensional Ising model. From the theoretical point of view, renormal- 
isation group calculations provide quite accurate results. Experiments carried 
out on crystalline mixtures of compounds lead to measurements as accurate as 
three digits on the values of critical exponents. Numerically, extensive Monte 
Carlo simulations then pretend to be of comparable accuracy. Life becomes 
complicated when details are compared between the three approaches. 
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1. Introduction 



Since many years the random Ising model has served as a paradigmatic system in 
which the influence of disorder may be studied through different techniques. Here we 
now mention three of them. The renormalization group (RG) approach, experimental 
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measurement and Monte Carlo simulation. 

RG calculations were considered quite early [1,2] and since then many groups 
were illuminated by skillful RG calculations. Among them, we would like to mention 
the work of Folk, Holovatch and Yavors'kii (for recent reviews, see Refs. [3,4]). On 
the experimental side, measurements on crystalline mixtures of difluoridc of different 
transition metals, e.g., magnetic FeF2 substituted with non-magnetic ZnF2, were 
performed in the same period over two decades, gaining in refinement and accuracy 
(see early results of Birgeneau et al. in Ref. [5], for a review sec, e.g., Ref. [3]). For 
the third aspect of simulations, technical progress made by computer manufacturers 
enabled more and more accurate simulations (which started for disordered systems 
in 3D, e.g., with Landau in Ref. [6]) and the study of disordered magnetic systems 
benefited from the development of parallel computing. Monte Carlo simulators thus 
competed in performance (for a review, see, e.g., Ref. [3]). 

To introduce the subject, we may imagine a conversation between three people 
a few decades ago, when the concept of universality was not as clearly stated as it is 
nowadays. Imagine a theoretician as the one who stands up with enthusiasm for this 
new concept. Let us call him Salviati. He has an interesting discussion with a good 
physicist, an experimentalist, let us say named Sagredo. It is usually considered that 
the concept of "numerical experiments" originates in the FPU problem, a numerical 
study of the thermalization of a chain of atoms, performed by Fermi, Pasta and 
Ulam at Los Alamos [7]. A third person is thus participating in our conversation. 
Having no more character free, let us call him Simplicio - the simulator! 

In a pub, beginning of the seventies: 

Sagr. Dear friends, I would like to report on recent experiments that I am conduct- 
ing at the lab. I produced many samples of difluoride of magnetic transition 
metals, substituting randomly different amounts of non-magnetic metal, and 
found very interesting results. 

Salv. It would be great to compare your results with recent theoretical predictions. 
Which quantity do you measure? 

Sagr. Critical temperature, correlation length, susceptibility,. . . Of course, the tran- 
sition temperature decreases when impurities are added, but what looks inter- 
esting is the neighbourhood of the transition. The singularities of some quanti- 
ties (susceptibility, correlation length) seem to be independent of the impurity 
concentration. 

Salv. This is a wonderful observation. It supports the universality assumption. 
You know, from recent RG theory, one expects that the free energy density 
has a singularity in the vicinity of the transition, and that this singularity is 
described by some critical exponents which are believed to be independent of 
the details of the system under consideration. 

Sagr. Do you mean that the presence of impurities is a detail? Experimentally it 
is not. It produces an observable decrease of transition temperature. 
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Salv. You are right, critical temperature is not universal and disorder is perhaps 
not a detail, but the precise amount of disorder probably does not matter, at 
least in some range. 

Sagr. But the singularities that I measured are different from those of a pure 
sample that I have also produced. A colleague of mine made similar experiments 
on Heisenberg-like samples and he did not notice any similar modification of 
the singularities due to the introduction of disorder. 

Salv. Probably that disorder is relevant in your case and not in his case. 

SiMPL. Maybe we could make Monte Carlo simulations, I have access to a computer 
and I have been told that it is not very difficult to produce simulations of a 
disordered Ising model on a cubic lattice. I only have to add quenched vacancies 
in the system and average over the disorder realizations. 

Salv. We may also write an effective (p'^ -theory for the diluted problem with a scalar 
field Ginzburg-Landau- Wilson Hamiltonian and calculate the critical exponents 
analytically. 

Sagr. / will measure many other quantities and see what is universal and what is 
not, and we will compare our results. 

And this is where the problems occur. In the comparison . . . 

We will give in the rest of the paper a short review of some recent progress in 
studies of the 3D disordered Ising model, emphasising the role of universality and 
its difficult emergence when trying to reconcile theoretical, experimental, and com- 
putational predictions. Reference will be made to seminal papers and to exhaustive 
reviews only. It is of course easy to wander a bit on arXiv and look around the 
names of Calabrcsc, Pelissetto and Vicari, Prudnikov, Shalaev or Sokolov, Folk or 
Holovatch. We apologise to those whose work is not directly mentioned in the short 
reference list, an indelicacy only due to our ignorance, our misunderstanding or our 
laziness - or all together. 



2. RG calculation of critical exponents 

Long distance properties of the Ising model near its second-order phase transition 
are described in field theory by an effective Ginzburg-Landau- Wilson Hamiltonian 



r) + f 



where ml is the bare coupling proportional to the deviation T — Tc from the critical 
point and </7(r) is a bare scalar field. Quenched randomness is introduced in such 
a model by considering that the adjunction of disorder results in a distribution 
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of local transition temperatures, so that a random temperature-like variable A is 
simply added to ttIq, 



H[if, A] 



(2) 



where A is drawn from, e.g., a Gaussian probability distribution of zero mean and 
dispersion o"^, V{A) = (27ro"^)~^ exp(— A^/2cr^). For a specific disorder realization 
[A], the partition function and the free energy read as Z[A] = J T)[ip]e^^^^'^''^^ and 
F[A] = — In Z[A]. Average over quenched disorder then requires to calculate 
quantities like F[A] = -P'^ J V[A]ln Z[A]V{A). This is performed through the 
introduction of n replicas of the model (labelled by a). Averaging over quenched 
disorder one ends up with an effective Hamiltonian with cubic anisotropy where the 
replicas are coupled through a new parameter vq 



replicas 



Tl Tl / 7X 

(x=\ (x=\ \a=l 




(3) 

Here the bare coupling uq, proportional to Mo, is positive and the bare coupling 
Wo, proportional to — o"^, is negative. The properties of the random Ising model are 
recovered while taking the limit n — > 0, InZ = lim„^o('Z^" ~ l)/''^- Under a change 
of length scale by a factor the field and couplings are renormalised according to 



.1/2 



0, mo 



rr 2 e-^U ^'-'V 



(4) 



where e — A 
parameters. 



D. The RG functions are defined by differentiation at fixed bare 



du 



9 In// 
dlnZ^ 



Pv{u,v) = 



dv 



d In ij, 



9 In// 



91nZ^2 
Q din II 



(5) 
(6) 



The skill of the theoretician is measured as his ability to compute these functions 
perturbatively, disentangling Feynman loops (they are known up to 6 loops), re- 
moving divergences which occur in the asymptotic limit by controlled rearrange- 
ment of the series for the vertex functions. Eventually expecting rehable results 
after complicated resummation procedures [8]. Fixed points are then solutions of 
(3u{u*,v*) = f3v{u*,v*) = 0, the stability of which is controlled by a stability matrix 
1^ with eigenvalues which besides the standard critical exponents also govern the 
corrections to scaling (exponent u). 

At that point, even Simplicio may read off the critical exponents! Consider for 
example the pair correlation function of bare fields {(p{0)(p{r)) . Under a change of 
length scale /(, it renormahses to Z^(/t)(0(O)0(r)). In the same manner, for another 
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dilatation parameter, /is, one has {(f{0)ip{sr)) — >• Z0(//s)(0(O)0(sr)). The ratio from 
this latter to the previous expression leads to 

(0(O)0(r)) imvir)) ■ ^ ^ 

This expression gives the algebraic decay of the two-point correlation function of 
renormalised fields (0(O)0(r)) ~ |r|-(-D-2+»7,A) in terms of the pair correlation func- 
tion of the bare fields which are described by mean-field theory (MFT), i.e., at 
the Gaussian fixed point (FP), {ip{0)ip{r)) ~ |r|-(-°-2) (t^mft = 0). The ratio 



din * 

Jfis^'i' evaluated at the new FP gives s~^'i' and leads to 



(0(O)0(r)) 

from which one reads off the value of the critical exponent at this FP: 

= ll- (9) 

Following the same argument, the scaling dimension 1/z/ of the (renormalised) tem- 
perature field w? is given at the random fixed point in terms of the MFT value, 
= 2, and one gets 

^ = 2-7;;- (10) 

From these two exponents, the others may be deduced by scaling arguments, de- 
scribing the leading singularities of the physical quantities, e.g., of the magnetic 
susceptibility: 

x{t)^T^\t\-\ 7 = 1/(2-77^). (11) 

In the non- asymptotic regime, the system approaches criticality in a more com- 
plex way and this is where corrections to scaling appear, 

x(t) ~ r±|T|-^(i + r'iirr'^ + r"±|T|''^'^ + ...), (12) 

where the scaling dimension uj corresponds to the negative of the leading irrelevant 
RG eigenvalue, u — — jyaj, as it is usually denoted, and the dots in Eq. (12) stand 
for higher order irrelevant corrections. Non solum, the critical exponents, sed etiam 
combinations of critical amplitudes and correction-to-scaling exponents are universal 
quantities. 

Also it is common practice, in order to describe the approach to criticality espe- 
cially in experiments and simulations, to introduce effective exponents through 

X(r)^r±|r|---, = (13) 

These effective exponents may be calculated theoretically from the flow equations, 
e.g., 

Vesil) = 7<^(^i(0,^(0), (14) 

= 2-Jrn^{u{l),v{l)). (15) 
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Figure 1. RG flows in the {u,v) plane. The points denoted "Gaussian FP", "Pure 
IM" and "RIM" arc respectively the Gaussian fixed point, the pure Ising model 
fixed point and the random Ising model fixed point. 

The variation of effective exponents depends on the RG flow in the parameter space 
as shown in Fig. 1. 

In the case of Heisenberg-like ferromagnets, the experimental observation of a 
maximum of the effective exponent 7eff found a theoretical explanation in terms of 

trajectories in the parameter space [9]. The same observation holds in the case of 
the random Ising model, but the critical exponents also change at the disorder fixed 
point in this latter case. 

3. Experiments 

Experiments on site-diluted three-dimensional Ising magnets are usually per- 
formed on uniaxial disordered anti-ferromagnets such as Fei_a,Zna,F2 or Mni_a,ZnxF2- 
The original aim was the study of the random-field behaviour when a uniform mag- 
netic field is applied to such a disordered system. However, when the samples are of 
high quality (low mosaicity, high chemical homogeneity), also the behaviour in zero 
external magnetic field is accessible {3D disordered Ising model universality class). 
Staggered susceptibility and correlation length are deduced from neutron scatter- 
ing experiments. The scattering intensity /(q) is the Fourier transform of the pair 
correlation function, where long-range fiuctuations produce an isotropic Lorentzian 
peak centred at the superstructure spot position qo with a peak intensity given by 
the susceptibility and a width determined by the inverse correlation length, while 
long-range order gives a background proportional to the order parameter squared: 




(16) 
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Fitting the Lorentzian at different temperatures eventually give access to the critical 
exponents, critical amplitudes, and possibly the correction to scaling. 

4. Monte Carlo simulations 

The majority of numerical studies of the disordered Ising model were concerned 
with site dilution. But we may also choose to model the disorder by bond dilution in 
order to compare these two kinds of disorder and to verify that they indeed lead to 
the same set of new critical exponents, as expected theoretically by universality. In 
our study we therefore considered the bond-diluted Ising model in three dimensions 
whose Hamiltonian with uncorrelated quenched random interactions can be written 
(in a Potts model normalisation) as 

-Pn = J2K^A.,'r,: (17) 

where the spins take the values (jj = ±1 and the sum goes over all nearest-neighbour 
pairs The coupling strengths are allowed to take two different values Kij — 

K = J/ksT and with probabilities p and 1 — p, respectively, 

nKrA = n piKrj) = n b^i^.. -k) + {i- p)6{K,,)i (is) 

c = 1 — p being the concentration of missing bonds, which play the role of the 
non-magnetic impurities. 

The phase diagram and the critical properties at a few selected dilutions were 
studied by large-scale Monte Carlo simulations on simple cubic lattices with V = 
spins (up to L = 96) and periodic boundary conditions in the three space direc- 
tions, using the Swendsen-Wang cluster algorithm for updating the spins. All phys- 
ical quantities are averaged over 2 000 - 5 000 disorder realisations, indicated by a 
bar (e.g., x for the susceptibihty) . Standard definitions were used, e.g., for a given 
disorder realisation, the magnetisation is defined according to m = where (. . .) 
stands for the thermal average and fi = {N^ — N^) / {N-^ + N^) with N-^^i counting the 
number of "up" and "down" spins. The susceptibility follows from the fiuctuation- 
dissipation relation, x = -^^((/^^) ~ (l/^l)^)- The phase diagram is obtained by 
locating the maxima of the average susceptibility xl (a diverging quantity in the 
thermodynamic limit) for increasing lattice sizes L as a function of the coupling 
strength K. 

As a function of the reduced temperature r = {Kc — K) {t < in the low- 
temperature (LT) phase and r > in the high-temperature (HT) phase) and the 
system size L, the susceptibility is expected to scale as: 

x{r,L)^\Trg±{L'/-'\r\), (19) 

where g± is a scaling function of the variable x — V-/^\t\ and the subscript ± 
stands for the HT/LT phases. Recalling (13) we can define a temperature dependent 
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Figure 2. Variation of the temperature dependent effective critical exponent 
7efr(|7"|) = — dlnx/dln |t| (in the low-temperature phase) as a function of the 
reduced temperature |t| (top) and V-/'^\t\ (bottom) for the bond-diluted Ising 
model with p = 0.7 and several lattice sizes L. The horizontal solid and dashed 
lines indicate the site-diluted and pure values of 7, respectively. 



effective critical exponent 7efr(|T|) = — dlnx/dln |r|, which should converge towards 
the asymptotic critical exponent 7 when L — 00 and |r| — 0. Our results for 
p = 0.7 are shown in Fig. 2. For the greatest sizes, the effective exponent 7efr(|T|) 
is stable around 1.34 when \t\ is not too small, i.e., when the finite-size effects are 
not too strong. The plot of 7cfr(|T|) vs. the rescaled variable L^^'^\t\ shows that the 
critical power-law behaviour holds in different temperature ranges for the different 
sizes studied. From the temperature behaviour of the susceptibility, we also have 
directly extracted the power-law exponent 7 using error weighted least-squares fits 
by choosing the temperature range that gives the smallest x^/d.o.f for several system 
sizes. The results are consistent with 7 ^ 1.34 — 1.36. 

From the previous expression of the susceptibility as a function of the reduced 
temperature and size, it is instructive to plot the scaling function g±{x). For finite 
size and |r| 7^ 0, the scaling functions may be Taylor expanded in powers of the 
inverse scaling variable = (L^^'^\t\)~^, x±{t, L) = \T\~"'[g±{oo) + x~^g'^{oo) + 
0(x~^)], where the amplitude 5'±(oo) is usually denoted by T±. Multiplying by L''^^" 
leads to 

X^L-^/" = x-^g±{x) = r±x-^ + 0{x-^-^). (20) 
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Figure 3. Log-log plot of the scaling function \L~'^I^ in the LT and HT phases 
(respectively top and bottom) against V^I^\t\ for p = 0.4, 0.55, and 0.7. The 
solid lines show the power-law behaviours with the exponent 7 ~ 1.34 which 
characterise the disordered fixed point. The data for the smallest values of L^/''|t|, 
which do not lyc on the master curve, illustrate the finite-size effects when the 
correlation length ^ is limited by the linear size L. 



The curves in the ordered and disordered phases, shown in Fig. 3, are obviously 
universal master curves whose slopes, in a log-log plot, give the critical exponent 
7 ~ 1.34. Indeed, when |t| — > but with L still larger than the correlation length 
^, one should recover the critical behaviour given by g±{x) — The critical 

amplitudes F-i- follow. 



5. Results and conclusions 

In conclusion, according to Sagredo's suggestion wc may compare the results 
deduced from the different techniques. We only concentrate here on the behaviour 
of the susceptibility, which already leads to partially conflicting results as can be 
seen by inspection of Table 1. How to conclude in favour of universahty? A possible 
answer would be the following. 
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In the same pub, a few days later, after inspection of the results: 

Sagr. Universality is still a good idea, but it is very difficult to produce high-quality 
samples where universality is clearly satisfied. 

SiMPL. For the simulations, it is so time consuming to increase the size of the 
system that it is not feasible at the moment. The problems here come essentially 
from the thermodynamic limit and the disorder average. 

Salv. In the equations, the sample is perfect, the disorder average is exact and the 
thermodynamic limit is automatically understood. So I believe that theoretical 
results are correct. 

SiMPL. AND Sagr. But what about RG calculations at 7 loop approximation? Will 
it come from St- Petersburg, from Roma, or from the Linz-Lviv axis? 



Table 1. Critical exponents and critical amplitude ratio of the susceptibility as 
measured with different techniques. 



Technique 


7 


r+/r_ 


UJ 


Ref. 


RG 




2.2 




[10] 




1.318 




0.39(4) 


[11,12)1 




1.330(17) 




0.25(10) 


[4]2 


Neutron scattering 


1.44(6) 


2.2 


0.5 


[5]^ 




1.31(3) 


2.8(2) 




[13]^ 




1.37(4) 


2.40(2) 




[14]5 


MC 


1.342(10) 




0.37 


[15]^ 




1.34(1) 


1.62(10) 


undetermined 


[16]^ 




1.342(7) 




undetermined 


[17]« 



1 4 loop approximation. 

^ 6 loop approximation, fixed dimension. 

3 Fei_^Zna;F2, x = 0.4, 0.5, |r| ~ lO'^. 

4 Feo.46Zno.54F2, 1.5 X 10-3 < |r| < W-\ 

^ Mno.75Zno.25F2, 4 X 10"" < |r| < 2 X 10-1. 
*5 site dilution, p = 0.4 to 0.8. 

bond dilution, p = 0.7. The correction to scaling is too small to be determined. 
^ site dilution, p = 0.8. The observed correction to scaling could be the next-to-leading. 
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